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Looking inside neutron stars: Microscopic calculations confront observations 
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While QCD appears not to be accurately solvable in the regime of interest for neutron star physics, 
microscopic calculations are feasible at both low and very high densities. In this work, we propose 
using the most realistic calculations in these two regimes of nuclear physics and perturbative QCD, 
and construct equations of state by matching the results requiring thermodynamic consistency. We 
find that the resulting equations of state — in contrast to several hadronic ones — are able to 
reproduce current observational data on neutron stars without any fine tuning, and allow stable 
hybrid stars with masses up to 2.1 M®. Using recent observations of star radii, we perform a 
maximum likelihood analysis to further constrain the equation of state, and in addition show that 
the effects of rotation on radii and masses should not be neglected in future precision studies. 



Introduction What are the thermodynamic proper- 
ties of matter found in the center of a pulsar? De- 
spite several decades of research in nuclear physics and a 
wealth of astrophysical data, this seemingly simple ques- 
tion has yet to be answered convincingly. It is well known 
that pulsars are made out of dense nuclear matter, and 
that the relevant interactions are described by Quantum 
Chromodynamics (QCD), but — unlike at high temper- 
ature and low density — accurately solving for the Equa- 
tion of State (EoS) at high densities has proven to be im- 
possible up to the present day. The situation is further 
complicated by the rich phase structure of cold nuclear 
matter, with possible phases containing e.g. kaon/pion 
condensation, hyperons and crystalline color supercon- 
ductors [2-|j|. For these reasons, a commonly used ap- 
proach in the study of high density nuclear matter is to 
give up on first principles calculations and resort to phe- 
nomenological models, such as those of the Nambu- Jona- 
Lasinio type (see e.g. Ref. [1] and references therein). 

In this article, we propose a different strategy, based on 
the fact that accurate microscopic calculations for bulk 
thermodynamic quantities can be performed for densities 
up to the nuclear saturation density n sat = 0.16 fm~ 3 in 
the hadronic phase 0, @ and at very high densities 
n 3> n sa t in the quark-gluon phase Q- It is known that 
for sufficiently high density, a deconfinement transition 
from the hadronic to the quark-gluon phase must take 
place [31j , and if one is interested in bulk thermodynamic 
quantities such as the EoS — in contrast to the detailed 
phase structure of nuclear matter — then it is plausible 
that quantitatively reliable results can be obtained by 
simply merging the high and low density results and re- 
quiring thermodynamic stability. Such an approach has 
indeed been found to be successful for nuclear matter at 
high temperature and low density, where a comparison 
to lattice QCD results is possible Q. 

In our work, we follow the general strategy outlined in 



Ref. [7|, generalizing its results to be applicable for neu- 
tron stars that rotate and/or contain a mixed phase of 
nuclear and quark matter. While attempts to constrain 
the nuclear matter EoS based on astrophysical data have 
recently been carried out [1, 0] , our work constitutes the 
first direct comparison of state-of-the-art theoretical pre- 
dictions with observations of rotating neutron stars. 

Methodology We begin by providing a brief account 
of how candidate EoSs for neutron stars can be obtained 
from microscopic calculations. Readers interested in the 
details of the setup are kindly asked to consult Ref. [7| . 

First, we must identify the most reliable and accurate 
results available for the hadronic and quark matter EoSs. 
The primary uncertainty in the hadronic sector comes 
from the composition of the matter, which could include 
nucleons Q , nucleons and hyperons Q , or nucleons with 
a kaon condensate @. Less significant uncertainties en- 
ter e.g. from neglecting the interactions between more 
than three nuclei, the form of the variational ansatz, the 
unknown form of the condensation potential, as well as 
from the hyperon interactions, which we fix to the values 
considered realistic by the authors of Refs. @, @. Our 
conclusion is that to obtain a conservative estimate for 
the uncertainty in the hadronicphase, it suffices to con- 
sider the three EoSs of Refs. as the limiting cases. 
For the very lowest densities, n < n sat /2, we use results 
from Refs. ^Ell- 
in the description of the quark-gluon phase, it has been 
customary to use the MIT bag model [li^ . which pro- 
vides an analytic EoS based on non-interacting massless 
quarks that is furthermore simple to parametrize. When- 
ever aiming for quantitative results, we, however, believe 
that it is not a fruitful approach to merge accurate (and 
complicated) calculations in the hadronic phase to the 
most simple-minded result in the quark-gluon phase, the 
popularity of which is to a large extent based on its ease 
of use. Thus, rather than using a model that altogether 
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FIG. 1: Various EoSs resulting from the matching of the 
hadronic and quark matter phases. Also shown are the pure 
hadronic EoSs for nucleons, nucleons+hyperons and nucle- 
ons+condensed kaons. Requiring hadron and quark matter 
phases to match in a thermodynamically stable manner sig- 
nificantly reduces the range of the allowed EoSs. 

neglects quark interactions, we employ the most accurate 
perturbative QCD calculation available, which includes 
interactions of quarks and gluons to three loop order as 
well as a running strange quark mass Q. The primary 
uncertainties that remain in the EoS are then the depen- 
dence of the result on the ratio of the renormalization 
scale to the baryon chemical potential A//ib (which nev- 
ertheless vanishes order by order in perturbation theory) 
and the integration constant B, which is to some extent 
analogous to the bag constant of the MIT bag model. 

Lesser uncertainties originate from the value of the 
strange quark mass m s , the scale parameter of the MS 
scheme A^jg as well as the gap parameter A, which en- 
ters the calculation should the quark-gluon matter be 
in a color superconducting phase. We have chosen to 
vary 3A//ib around its canonical value 2 by a factor 
of two, and let B take any value that is not in con- 
flict with the positivity of the energy density or ther- 
modynamic stability (c/. Ref. [7|). The parameter m s 
is fixed using the most recent lattice results indicating 
m s (A = 2GeV) = 92.2 ± 1.3 MeV 0, while A^ and 
A are let to vary in their expected ranges 346 MeV 
< A^s < 412 MeV Q and < A < 100 MeV 0. We 
stress that with A = 100 MeV, we probe the maximal 
effect that color superconductivity can have on the EoS. 

While both the hadronic and quark matter calculations 
come with considerable uncertainty bands, it turns out 
that when requiring thermodynamically stable matching 
between the two phases, the uncertainty is significantly 
reduced (see also Fig. Q}. We find that this is true 
both for homogeneous phases and when allowing a mixed 
phase of quarks and hadrons, the latter scenario being 
viable if the microscopic surface tension of QCD is low 
enough 
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In both cases, we were unable to match 
the quark phase onto hadronic EoSs where hyperons are 
present or kaons have started to condense, because the 
pressure and energy density of these EoSs exceed those of 



Name 


Spin (Hz) 


Mass/M 


Remarks 


J0024-7204H 


312 


1.41 ±0.08 


in 47 Tuc, [16] 


J0437-4715 


174 


1.76 ±0.20 


[15, 2^ 


J0514-4002A 


126 


< 1.52 


in NGC1851, [22J 


J0751+1807 


288 


1.26 ± 0.14 


[22] 


J1012+5307 


190 


1.64 ±0.22 


[16] 


J1713+0747 


219 


1.53 ±0.08 


[15, 17] 


4U1608-52 


619 


1.70 ± 0.40 


LMXB, [8, 23] 


J1748- 24461 


105 


1.85 ±0.05 


in Ter 5, [22] 


SAXJ1808.4-3658 


401 


1.40 ± 0.20 


LMXB, [26] 


J1824-2452C 


240 


< 1.37 


in M28, [22] 


B1855+09 


187 


1.58 ±0.13 


[16] 


J1903+0327 


465 


1.67 ±0.01 


[15, 25] 


J1909-3744 


339 


1.44 ±0.02 


[15, 17] 


J1911-5958A 


306 


1.40 ± 0.16 


in NGC 6752, [22] 


J2019+2425 


254 


< 1.51 


[16] 



TABLE I: Masses of neutron stars with millisecond periods; 
see Refs. @, [H4T3, dJii] and references therein. 

quark matter at the same chemical potentials. Consider- 
ing that at extremely high densities, nuclear matter must 
be in the quark-gluon phase, we take this as an indication 
that neither hyperons nor kaon condensation are realized 
in cold nuclear matter. As will be demonstrated below, 
the same conclusion can be independently reached based 
on neutron star observations. 

Results Neutron star observations have progressed 
significantly since the first discovery of a neutron star 
more than 40 years ago. The spin frequencies of more 
than 1800 objects have been catalogued to date [l5j|, the 
masses of approximately 70 neutron stars are known with 
some precision [3, [ItJ , and very recently, observational 
constraints on neutron star radii have become available 

am 



321 ] . In the following, we will use these ob- 
servations in an attempt to constrain the cold nuclear 
matter EoS. Neutron stars with millisecond periods, for 
which accurate mass determinations exist, will turn out 
to be particularly useful in this process, and to this end 
we have collected all 15 presently known objects in Tab.fl] 
While it is customary to use only mass-radius data to 
constrain the cold nuclear matter EoS, we emphasize that 
observations on neutron star masses and frequencies offer 
an interesting alternative because of the relatively accu- 
rate data available for many such systems (see Tab. [TJ) . 
This is illustrated in Fig. [21 where the minimal and max- 
imal neutron star masses allowed by the solutions to the 
general relativistic field equations for uniformly rotating 
stars are shown for a class of different EoSs. From here, 
it is seen that the EoSs containing hyperons and con- 
densed kaons imply maximal masses for the stars that 
are significantly below the masses of several observed ob- 
jects. This strengthens the above finding that hyperons 
and condensed kaons are most likely not realized in cold 
nuclear matter. 

For EoSs consistent with observational data on mil- 
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FIG. 2: Theoretical upper and lower limits for the masses 
of rotating neutron stars plotted against observational data. 
If the observed stars are approximately uniformly rotating, 
EoSs with hyperons and condensed kaons can be ruled out. 



Name 



Mass/Me Radius/km 



4U1608-52 

4U1742-307 

EX01745-248 



1.70 ±0.43 11.5 ±2.0 
1.63 ±0.58 14.4 ±1.0 
1.55 ±0.40 10.0 ±1.8 



SAXJ1808.4-3658 1.40 ±0.20 10.0 ± 4.0 
4U1820-30 1.73 ±0.33 11.5 ±1.5 

RXJ1856-375 1.70 ±0.30 11.5 ±1.2 

TABLE II: Known neutron star masses and radii [^.IsL fTsL l20j| . 

lisecond period stars, we point out that stability with re- 
spect to radial oscillations as well as mass-shedding im- 
pose further limits on the masses, and that this range 
is rapidly shrinking for the fastest spinning objects. The 
maximum masses tend to increase markedly with the spin 
frequency /, for some EoSs as much as 20 percent in com- 
parison with / = 0. On the other hand, using the the- 
oretical minimum mass for the iron core of a progenitor 
star as input [27J , current supernova simulations suggest 
that newborn neutron stars have a minimum (gravita- 
tional) mass of ~ I.IMq 28[ (also indicated in Fig. [5]). 
As a consequence of the frequency dependence of these 
lower and upper limits, we expect stars with higher spin 
rates to have higher masses on average. Improvements in 
our understanding of core-collapse supernovae together 
with future observations could make it possible to even- 
tually accentuate the conclusions to be drawn from Fig. [2] 
For those EoSs that have not been ruled out above, 
it is possible to continue the comparison with available 
mass-radius data. At present, neutron star radii have 
been determined with reasonable accuracy for altogether 
6 stars (Tab. HQ, so that for two of them, 4U1608-52 and 
SAXJ1808. 4-3658, information on three different quan- 
tities — mass, radius and spin — is available. Neutron 
star radii are customarily determined assuming the stars 
to be non-rotating. If we follow this assumption and per- 
form a maximum likelihood analysis of the six objects in 
Tab. UH we find that EoSs generating stars with a large 
quark core and only a small hadronic crust are favored. 
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FIG. 3: Star mass-radius curves for hybrid EoSs compared 
to observations (both assuming negligible spin frequencies). 
The pure phases region is split into two parts according to 
the 'Case F and 'Case IF matching criteria of Ref. 0. 

The likelihood is measured in terms of a merit function 
X 2 = Eli Mm[(M-M 2 )>L± {R-Ri) 2 /a 2 Ri \, where 
Mi, Ri are the masses and radii from Tab. [TTl omi and 
am are the respective standard deviations, and the min- 
imization is performed along the MR curve generated 
by the EoS. We observe that the large central value and 
small error bars for the radius of 4U1724-307 have the ef- 
fect of disfavoring EoSs with small radii {e.g. those with 
mixed phases) and favoring ones where the phase transi- 
tion takes place at moderately small densities. 

Overall, all of the hybrid EoSs that we obtain by 
matching the hadron and quark phases result in mass- 
radius curves that are reasonably close to observations, 
as can be seen in Fig. [21 For slowly rotating hybrid stars, 
masses up to 2.1 Mq are possible, but these stars con- 
tain only a tiny quark core at the center of the star. 
When considering a mixed phase of hadrons and quarks, 
the maximal masses encountered are 1.93 Mq, and for 
hybrid stars with only a small hadronic crust, we find 
masses less than 1.9 Mq. In general, for hybrid stars 
with masses of 1.5 Mq, the radii range from 10.8 to 13.5 
km, with 13.1 km corresponding to the smallest \ 2 value 
when attempting to fit to observations (c/. Fig. [3]). 

Finally, we note that the assumption of non-rotation in 
the analysis of neutron star radii introduces systematic 
errors, as many of the stars are in fact rapidly rotating. 
In Fig. |4l we demonstrate how the polar and equatorial 
radii of a star with a constant baryonic mass start to 
differ strongly at large frequencies. When aiming at sub- 
km uncertainties in radii observations and theory analysis 
in the future, the effect of rotation will thus have to be 
taken into account for frequencies exceeding 600 Hz. 

Conclusions We have constructed possible equations 
of state for neutron star interiors by matching state-of- 
the-art microscopic calculations at low and high densities 
(hadronic and quark matter phases, respectively). Our 
findings show that these EoSs lead to neutron star masses 
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A/MeV 3A/(^isA 5T g) x GeV M match /MeV mixed phase \ 2 



100 


9.26 


955 


no 


4.42 





11.56 


960 


no 


5.93 


100 


11.56 


1430 


yes 


9.07 





8.67 


1617 


yes 


9.16 


100 


10.58 


1222 


no 


9.70 





11.56 


1492 


no 


10.86 





11.56 


1492 


yes 


10.87 


100 


11.56 


1257 


yes 


14.48 



TABLE III: \ 2 values for selected hybrid EoSs with quark 
matter (parameters A, 3A//i_g and Ajjg) matched to nucleons 
at ub = Mmatch- The set consists of the best and worst fits 
for A = 0, 0.1 GeV and with and without mixed phases. For 
the purely nucleonic equation of state, we obtain \ 2 = 9.16. 




FIG. 4: Observational data for the radii and frequencies of the 
stars SAXJ1808.4-3658 and 4U1608-52. Using the minimum 
X 2 hybrid EoS (A = 0.1 GeV, 3A/( A i s A S s) = 11.56 GeV" 1 , 
Mmatch = 956 MeV, no mixed phase) and assuming a star mass 
of M/Mq ~ 1.5, results for both the polar and equatorial radii 
as functions of frequency are shown for comparison. 

and radii that are consistent with present observations 
without having to fine tune any parameters. Using cur- 
rent observational data, we are able to rule out hadronic 
equations of state that contain either hyperons or as- 
sume kaon condensation. A maximum likelihood anal- 
ysis seems to favor neutron stars consisting of a large 
quark matter core with only a thin hadronic crust, but 
— given the scarcity of objects, for which sufficient data 
is available — this conclusion must be considered tenta- 
tive. Finally, we have pointed out how to further con- 
strain the EoS once additional observational information 
on neutron star masses, frequencies and radii becomes 
available, in particular for rapidly rotating systems. 

All EoSs used in this article are available online [3fj| . 
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